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We report on preliminary results of a high statistics quenched lattice QCD calculation of nucleon matrix 
elements within the Symanzik improvement programme. Using the recently determined renormalisation 
constants from the Alpha Collaboration we present a fully non-pertubative calculation of the forward 
nucleon axial matrix element with 0{a) lattice artifacts completely removed. Runs are made nt (3 — 6.0 
and /3 = 6.2, in an attempt to check scaling and 0{a^) effects. We shall also briefly describe results for 
(x), the matrix element of a higher derivative operator. 



1. INTRODUCTION 

In this talk we shall describe results for nu- 
cleon matrix elements: (A^|0-^|A^) (at p — 0) 
using 0(a) Symanzik improved fermions for 



the vector current, 
warm-up exercise) 



(iV|V^|iV) (as a 



• the axial current, {N, s\Aj^\N, s) — 
Ags^, where Aq is the fraction of the 
nucleon spin carried by the quark q 

• (x) , the fraction of the nucleon mo- 
mentum carried by quark q 

We have worked in the quenched approxima- 
tion and generated 0(500) configurations at 
/3 = 6.0 on a 16^ x 32 lattice and 0(150) con- 
figurations at = 6.2 on a 24'' x 48 lattice. 
In both cases we have used three k values to 
enable us to extrapolate to the chiral limit. 
For the hadron spectrum see Q . The method 
to determine the matrix elements is standard. 



see eg ; we only note that we are computing 
just the quark line connected term. 

2. IMPROVED OPERATORS 

Symanzik improvement is a systematic im- 
provement of the action and operators to 
0(a") (here 0(a^)) by adding a basis of irrele- 
vant operators to completely remove 0(a"~^) 
effects. Restricting improvement to on-shell 
matrix elements means that the equations of 
motion (EOM) can be used to reduce the set 
of operators. For the action we only need one 
additional operator - the clover term, ^ , with 
known coefficient Cswig^), 0- We write for 
the improved axial, vector currents: 



V. = 



Vfj^ - civaq Dp.q + 

-2m'j^{x), with 



*Talk given by R. Horsley at Lat97, Edinburgh, U.K. 



while {N\%,-^j:^'rjN) 
-C2aq D{f,Du} q 
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and Dfj,= |(-Dp ^ Dfj.)- Derivative opera- 
tors which do not contribute to the forward 
matrix element have been dropped. Rcnor- 
mahsed operators are then given by — 
Zo{l + boamq)0. Note that these 0{a) op- 
erator sets with coefficients bo, {c/} are over- 
complete - using the EOM they can always 
be reduced by one. 

3. MATRIX ELEMENTS 



For the vector current we have from cur- 
rent conservation (iV| V4'^''^| A'^) = x^^\ with 
X*") = 2 and x^'*' = 1- Using the EOM 
we have a one parameter degree of freedom, 
which we can take to be Cw- From linear 
fits (in amq = I/kJ) to {N\Vi\N), 

—a{N\q Di q\N) we can find Zy from the con- 
stant terms and from the gradients hv- (Lin- 
ear fits seem to be adequate within our quark 
mass range aniq < 0.1.) The result for Zy, by 
is shown in Fig. ^ for /3 = 6.0. The Alpha 
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Fi gure 1. Results for Zy and by against ciy for 
/3 = 6.0. The error bars are given by bands of dotted 
lines. Also plotted is the one-loop perturbative result 
for Zy and the Alpha Collaboration results. 

Collaboration has non-perturbatively found 
Zv{civ = 0) and 6y(ciy = 0), [|. We see 
from the figure that there is very good agree- 
ment between the different determinations for 
Zv{civ = 0) and although for bv{civ = 0) 
there is some discrepancy we note that this is 
not unexpected as by has 0{a) effects. Also 



shown are two points (stars) from a method 
described in Again there is good agree- 
ment (this time at a non-zero value of cw)- 

For the axial current, again the Alpha 
Collaboration has non-perturbatively found 
■^a(cia = 0), H). This is sufficient for us 
as we are only interested in the results in 
the chiral limit. Performing these extrapola- 
tions for /3 — 6.0, 6.2 gives the scaling plot. 
Fig. H To attempt a comparison with Wil- 
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Figure 2. A scaling plot of Auv — Arf„ = against 
, using the string tension to set the scale. The 
Symanzik improved data at /3 = 6.0, 6.2 is shown 
with filled squares, the Wilson data, open square, is 
at f3 = 6.0. The physical value of gj^ = 1.26. 

son data, we have also plotted the result at 
f3 = 6.0, 1^ using the non-perturbative Za 
as found in ||^ of 0.78. This is close to the 
value (0.75) given from linearly interpolating 
the non-perturbative Ward Identity results in 
to /3 = 6.0, which would indicate that 0(a) 
effects between these two methods are small. 
In Fig. ^ we expect to have to make a linear 
extrapolation in for the improved case. 

For (x) there is no non-perturbative de- 
termination of the renormalisation constants 
yet available. First order perturbation theory 
gives, however, numbers rather close to 1, 
and indeed using tadpole improved (TI) per- 
turbation theory, we see that with one deriva- 
tive operators the (large) tadpole terms can- 
cel. This might indicate that perturbation 
theory is reasonably correct. From the defi- 
nition of the improved operator we see that 
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we have parameters Bq, ci, C2. To remove 
0(a) effects completely we must only have 
a residual one parameter degree of freedom 
(from the EOM), so that, eg, C2 = C2(ci) and 
bo = &o(ci). These functions are unknown at 
present. At tree level C2 = ci, bo — 1 ~ ci, 
which we shall use here. In Fig. |^ we compare 
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Figure 3. A comparison between first order per- 
turbation tlieory and TI perturbation theory, using 
e boosted coupUng ci'^'^{^a) given in table I of 
Also used in the perturbative coefficient IS Csw — 

one loop perturbative results, and TI results 
for (a;)(") in the chiral limit. We know that 
upon using the true coefficients then the re- 
sult must be independent of ci. TI appears 
to achieve this somewhat better than 1-loop 
perturbation expansion, so we shall use this 
in our scaling plot shown in Fig. ^. This 
operator has an anomalous dimension, so we 
must scale the results to the same /i; this 
has been performed with the scaling formula 

. Scal- 
6.0 where 



ing the (3 = 6.2 result to (3 
jjL ~ 1.95(GeV)~^ gives a 4% increase. We 
compare with the phenomenological value, 
Our results seem rather constant (in a^), sug- 
gesting that even in the continuum limit the 
lattice value is too high. We can only specu- 
late on the discrepancy: most likely this is due 
to a lattice problem - quenching, chiral ex- 
trapolation, Z not accurately enough known 
or perhaps there is a phenomenological prob- 
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Figure 4. Scaling plot of non-singlet {x)]^' — (x) 
against , using ci = 0. Same notation as in Fig. ti. 

lem - the fits are made from low so there 
might be higher twist contributions, pO|. 
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